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Abstract: We report observations of the ro-vibronic 
a 1 A — > X 3 £~ transition of SO at 1.707 (im on Io. 
These data were taken while Io was eclipsed by Jupiter, 
on four nights between July 2000 and March 2003. We 
analyze these results in conjunction with a previously 
published night to investigate the temporal behavior of 
these emissions. The observations were all conducted 
using the near-infrared spectrometer NIRSPEC on 
the W.M. Keck II telescope. The integrated emit- 
ted intensity for this band varies from 0.8 x 10 27 to 
2.4 x 10 27 photons/sec, with a possible link to varia- 
tions in Loki's infrared brightness. The band-shapes 
imply rotational temperatures of 550-1000K for the 
emitting gas, lending further evidence to a volcanic 
origin for sulfur monoxide. An attempt to detect the 
B 1 !] — > X 3 £~ transition of SO at 0.97 (im was un- 
successful; simultaneous detection with the 1.707 /jm 
band would permit determination of the SO column 
abundance. 



1. Introduction 

Io, the most volcanic of Jupiter's moons, has been 
scrutinized from both ground-based telescopes and 
spacecraft fiybys for many years. In 2003, Galileo 
completed its extended mission to the Jovian system 
which provided a wealth of information, and an insight 
into the nature of Io's atmosphere and surface. Like 
most good science, the data raised many new ques- 
tions, as it answered old ones. This volatile moon 
provides us with an intriguing testbed for examining 
the effects of tidal stresses on massive bodies, while 



ever-changing atmosphere (jLellouch et alj 11996). It 
remains our most easily observed example of active 
volcanism in a low pressure environment, and of a dy- 
namically supported atmosphere; an analog which may 
be particularly relevant to Martian and Lunar history. 

The first detection of Io's atmosphere was made by 
the Pioneer 10 fl y by which detected an ionosphere near 



the terminator (jKliore et al.l 119741 ) . The first atmo- 
spheric species, sulf ur dioxide, was id entified in Voy- 
ager IRIS spectra bv lPearl et all (| 19791 ). and confirmed 



at m illimeter and UV wav e length s by iLellouch et al 



(|l990h and iBallester et"aH (|l994h respectively. The 



observations and subsequent analysis posed two ma- 
jor questions about the atmosphere which are yet to 
be fully resolved. What are the relative contribu- 
tions to the atmosphere from volcanic eruptions, sub- 
limation, and sputtering? And how localized is the 
atmosphere? It is widely agreed that all the above 
mechanisms must play a part although to what degree 
is not yet well understood. The presence of spatia l 
variations was confirmed by iMcGrath et al.l (l2000h . 
where they used spatially resolved HST observations 
to confirm an atmospheric density enhancement over 
Pele, a volcanic hotspot, compared to a quiescent loca- 
tion, "T3". This local atmospheric enhancement near 
volcanic plumes has sin ce been modeled in detail for a 
variety of plume types (jZhang et al. I l2003l) . Currently 
the exact properties and dominant processes of Io's 
global a tmosphere are still debated. Numerical mod- 



els Ce.g. IStrobel et a"Dll994L IWong and Johnson! \l99^L 



Wong and 




i ts co nt inuous volcanic activity (e.g . iMarchis et al 
( 2005 ). iRathbun and Spencer ( 2006^ 1) results in an 



provide an estimate of the 



vertical thermal and chemical profile; however, these 
models depend on the local surface temperature and 
pressure, which are not well constrained. 

While the SO2 component of the atmosphere has 
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been studied in depth, identification of other gases in 
the atmosphere has been slow. The initial detection o f 

(1996) 



SO came in the millimeter when Lellouch et al 



identified it's rotational lines at 220 and 138GHz. The 
first detection of ro-vibronic transitions of the sulfur 
monoxid e radical in the infrar ed at 1.707/im was re- 
ported in lde Pater et al.l ( 2002t ) (henceforth referred to 
as DP02). 

The importa nce of quantifying the l e vels o f SO was 
emphasized by ISummers and Strobell ( 19961 ). They 
suggested that SO's vapor pressure may be high 
enough that the gas may not condense ou t on Io's 
night side, which supported the predictions of IngersoU 
(1989). ..t could, along with O2, be an important buffer 



gas in the atmosphere, allowing nanobars pressures to 
exist following the collapse of the detected but con- 
densable SO2 atmosphere. However SO is expected 
to be rapidly removed from the atmo sphere through 
chemical reactions on the surface (Lellouch et al 



1996), and hence detections of this gas while Io is 
in eclipse imply that it must be constantly replen- 
ished. Volcanic production of SO via the quenching of 



the g ases in high pressure vents (jZolotov and Feglev 
Il998l i was shown to be able to support the observed 
ratios of 3-10% by volume of SO/SO2. Thus we should 
expect localized enhancements of SO abundance near 
such hotspots. 

The SO ro-vibronic infrared band was detected in 
emission in Io's near IR spectrum (1.707/im), while the 
satellite was in Jupiter's shadow. DP02 fit a rotational 
tempera ture of 1000K to the band using lab spec- 
tra from Setzer et al.l (Il999|) along w i th lin e strengths 
from Bellarv and Balasubramanianl (1987), support- 
ing Zolotov and Feglev! ( 1998h 's theory that the most 
likely source of this gas was a high temperature vent. 
This current paper, a follow up to DP02, identifies 
this band on additional nights, supports their rota- 
tional temperature estimates and presents evidence of 
variability in the band emission. 

2. Observations 

We present observations of Io taken while the satel- 
lite was in Jupiter's shadow on four nights, from 2000 
to 2003, and compare these to the 1999 data pre- 
sented in de Pater et al 2002. The data were taken 
using the near infrared spectrometer NIRSPEC on the 
W.M.Keck II telescope at Mauna Kea, Hawaii, which 
is equipped with a 1024 x 1024 InSb ALADDIN de- 
tector. The data presented here used the NIRSPEC- 
6 (N6: 1. 558-2. 315/im) filter in low resolution echelle 
spectroscopy (R=2000) mode. The spectral region of 
the N6 filter analyzed for this paper was 1.6 - 2.0 /im. 
The total exposure time in each case was a few min- 



utes, which led to the detection of the SO band at 
1.707/im on all of the nights, although, as we will 
show later, at various strengths. In all cases Io was 
observed during the entire time it was in Jupiter's 
shadow, which was close to an hour in duration on each 
night. The typical viewing geometry of Io is shown in 
Fig. [T] , which shows some of the hotspots visible on 
Io going into and emerging from eclipse (2003 data). 

A list of observations and settings is given in Table 
[TJ We used a slit size of 0.76" x 42", the widest avail- 
able, with 0.144" pixels. Io subtends ~1" and hence 
is larger than the slit. We also used a slit camera 
(SCAM) which has a 256 x 256 pixel HgCdTe PIC- 
NIC array of 0.18" pixels, which gives a field of view 
of 46" x 46". These images were used to confirm Io's 
position on the slit. One of the data sets was taken 
with the adaptive optics system, on 1 2th Nov 2002. 



These data are discussed in detail by Ide Pater et al 
(|2007l ). In this paper we report the disk integrated 



total flux density of these measurements. 

All data were reduced normally: flat fielded, bad 
pixels removed, calibrated on neon and argon lines, 
and fine-tuned with telluric OH. The resulting cali- 
bration is correct to < lA. To remove telluric emis- 
sion lines we subtracted 'nodded' pairs of spectra 
which were taken in quick succession. We then av- 
eraged 15 pixels over Io (approximately twice the full 
width at half maximum) and subtracted the sky back- 
ground, to obtain the highest signal-to-noise spec- 
trum possible. We corrected for telluric absorption 
lines by dividing by either a normalized G or A star 
spectrum. The spectra were then multiplied by an 
appropriate normalized stellar model to correct for 
the known stellar features. G stars are moderately 
flat in this region although they have many small 
features. Alternatively, A stars are smoother with 
only a few well known pronounced features which can 
be easily modeled. The models were also used for 
flux calibration and are discussed further below. All 
spectra were rectified and partially processed using 
routines from the REDSPEC pipeline developed for 
NIRSPEC reductions by the UCLA infrared lab, (see 



http : / / www2 . keck . hawaii . edu/inst /nir spec/ redspec / ) 
Further analysis was performed using IDL routines 
written by the authors. 

The flux calibration of the spectra proved to be a 
challenging problem, due largely to the fact that the 
largest slit width for NIRSPEC is 0.76" compared to 
Io's angular size of ~1". In DP02, this problem was 
tackled by calibrating Io's flux density in the SCAM 
(slit camera) images, by observing the G 4 infrared 



standard star SJ9182 (Per sson et al. Il998h . Io's con 



tinuum level in the reduced spectra was then scaled to 
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Table 1 

Observations of Io in eclipse taken between 1999-2003 



Date(UT) 


Time(UT) 


Object 


Int Time(s) 


SEL 1 (deg) 


Ang Size(arcsec) 


Airmass 


9/24/1999 


11:02 


Io in eclipse a 


240 


348 


1.236 


1.065 


9/24/1999 


09:18 


SJ9182 b 


60 






1.384 


7/19/2000 


14:18 


Io in eclipse 


150 


347 


0.91 


1.95 


7/19/2000 


12:55 


HD201941 


75 






1.1 


11/12/2002 


14:21 


Io in eclipse c 


90 


346 


0.98 


1.13 


11/12/2002 


12:53 


HIP44336 


180 






1.38 


12/21/2002 


12:44 


Io in eclipse 


360 


348 


1.098 


1.03 


12/21/2002 


13:14 


HIP 45874 


40 






1.01 


3/25/2003 


08:52 


Io in eclipse 


180 


13 


1.071 


1.13 


3/25/2003 


08:40 


SAO 98024 


120 






1.1 



1 Io's Geodetic Sub Earth Longitude 
published in DP02 



TR standard iPersson et al.1 (|1998l ) 



= Disk averaged AO data, see Ide Pater et al.l ((2007) for full details 



match the SCAM flux level, which is measured in the 
same filter bandpass and along the same optical path. 
The 2002 AO data from DP06 were also calibrated in 
this manner. The lack of suitable SCAM images taken 
on the other nights in 2000-2003 prevented a similar 
method being used. 

In this paper we directly calibrated the stellar spec- 
tra using 2MASS. We used stellar models created by 
Kurucz (http://kurucz.harvard.edu/stars), which we 
reddened using the Fitzpatrick 1999 extinction curve, 
and flux calibrated to match the 2MASS magnitudes. 
The telluric calibration stars, despite not being pub- 
lished IR standards, all have been measured by the 
2MASS survey, and all have the highest 'AAA' photo- 
metric stability rating. The photometric errors quoted 
for these stars are all < 5%, which is adequate as this 
is dwarfed by the uncertainties due to centering of the 
object on the slit and variations in atmospheric condi- 
tions. 

For the non-AO Io data we corrected for the over- 
filling of the slit by convolving a model disk of Io, 
with a model gaussian PSF with a width equal to 
the average seeing on each night. Io's angular size, 
which varies with distance to earth, ranged from 0.91" 
to 1.24" on the various nights. The average seeing 
was 0.50" ,0.53" ,0.65" and 0.70" for the 1999, 2000, 
2002 (Dec) and 2003 nights respectively. We could 
then estimate the percentage of light which entered 
the slit (69-80%) and adjust the photometry appropri- 



ately. The photometric errors after this correction are 
dominated by the variation of our calibration stan- 
dards due to the atmosphere and positioning of the 
targets on the slit. We estimate the uncertainty to be 
~ 25-30%. Obviously this is not ideal, however it is the 
best achievable given the data set, and is sufficient for 
our purposes. To ensure consistency we reanalyzed the 
1999 data, using the flux calibration method described 
above. In the original paper, DP02, the total photon 
flux of the integrated SO band was found to be 2.3 x 
10 27 photons/sec, compared to 2.4 x 10 27 photons/sec 
using the above method. The September 1999 data 
presented below has been re-calibrated using this slit- 
loss correction method and are shown for comparison 
with DP02. 

Finally, as the distance from Earth to Io varies 
from 4-6AU, for simple comparison of the fluxes from 
Io, all nights have been adjusted to the Earth-Io dis- 
tance(4.08AU) of the 1999 data from DP02. It should 



be no ted that the data presented in Ide Pater et al 



()2007l ) have not been adjusted to this distance hence 
the apparent discrepancy between the scales of the 
November 2002 AO disk integrated data in these two 
papers. The results are later converted to total emit- 
ted intensity (distance independent, see Section 13. 3p 
for ease of comparison with any future observations. 
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3. Analysis 



Io 

11:02 UTC 
24th Sep 1999 




Io 

8:52 UTC 
25th Mar 2003 




Fig. 1— Observed face of Io during 1999 and 2003 
eclipses. The dashed lines indicate degrees longitude 
and latitude. Some of the most prominent hotspots 
and the sub-earth point (SEP) are indicated. The 2000 
and 2002 eclipses present nearly identical faces as the 
1999 data, as these are all eclipses in ingress, whereas 
the 2003 eclipse is in egress. 



3.1. Continuum Analysis 

To analyze the emission band we subtracted a single 
temperature blackbody fit to the continuum of the sur- 
rounding spectra from 1.63-1.67^tm and 1.75-1.79/xto, 
wavelength regions least affected by telluric lines and 
Io's SO emission band. In reality of course this black- 
body emission is the superposition of many different 
temperature hotspots, any of which can vary over mul- 
tiple observations. There is no reason to believe a di- 
rect correlation between total blackbody emission and 
SO band strength exists, as the S O may be a product 
of on ly some types of eruptions (jZolotov and Feglev 
19981) . 

Measurements of Io's continuum are complicated 
by the fact that by their nature all eclipse measure- 
ments of Io are taken when the satellite is very close 
to Jupiter (20" -40") and thus Io's spectra can be con- 
taminated to varying degrees by Jupiter's glare. We 
therefore compare our results with a theoretical esti- 
mate of the variability of the blackbody spectrum in 
the wavelength range of the N6 filter. To do this we 
used the temperature s and surface areas calculated in 
de Pater et al.l (|2004T ) as the basis for a 19 hot spot 
model of Io. Although these values were calculated in 
the 2.5-3/Ltm range , where we are not as sensitive to 
high temperature hotspots, they are a sufficient start- 
ing point for this estimate. We used a 10 5 iteration 
monte carlo simulation where the temperatures var- 
ied with gaussian fluctuations, with a la value of 20% 
of the starting temperature. For each distribution of 
temperatures the photon flux can simply be estimated 
from: 



N 



A 



$ = 2^—B i {\,A i ,T i )cos( 



0cos(6»i)4- 



Where Bi is the emitted flux of the i hotspot 
which depends on th e temperatures and areas from 



de Pater et alJ (|2004l ). Bi and fa are the longitude and 
latitude of the i th hotspot w.r.t. the central meridian 
longitude and sub-earth latitude, and di is the geocen- 
tric distance to Io (4.08AU). The cosine factors include 
the effects of foreshorte ning for each ho t spot using the 
positions determined in lde Pater et al. (2004) and the 
typical sub-earth point for eclipse observations of Io in 
ingress (~347 degrees longitude). Fig. [2]shows the 1 a 
upper and lower bounds of the flux at each wavelength 
for 10 5 iterations. Over-plotted are the fitted continua 
from the data taken in 1999-2003. The observed flux 
densities all fall between the 1 a limits of the Monte 
Carlo model, as expected. 
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Fig. 2. — Above shows the 1 sigma limits of a Monte 
Carlo model showing expected flux with hotspot tem- 
peratures varying by 20%, over-plotted are the fitted 
background continua from each night. 



We also modeled a system in which all of the 
hotspots are stable except for one. We examined 
the effect of an outburst by varying Loki's temper- 
ature by up to 50%, and saw a very similar pro- 
file to that shown in fig [5J This implies that a sin- 
gle, larg e hotspot can domi n ate the thermal back- 
ground (iMarchis et alj l2002l iMacintosh et alj 120031 
de Pater et al.1 120071 ) . adding more uncertainty to the 
background measurements. We also investigated the 
effects of varying the hotspot size by 10% at constant 
temperature and varying the sub-earth point over the 
hour duration of Io's eclipse and found these to be of 
minor influence compared to the variations caused by 
temperature fluctuations. 



3.2. Rotational Temperature 

As reported in DP02, the emission band detected at 
1.707^m is the SO forbidden a 1 A — > A 3 S~ ro-vibronic 
transition. The high resolution sp e ctrum of this tran- 
sition at 322K from ISetzer et al. ( 19991 ) is shown in 
Fig. [UJa). In panel (b) we show the same band after 
convolution with a gaussian of full width at half maxi- 
mum of 8.5A, which very closely approximates the line 
spread function for NIRSPEC. 

To investigate the emission band's dependence on 
rotational tempera ture, we use DP02's model, with 
wavenumbers from ISetzer et all ( 19991 ) and the rota- 
tional line-strength factors for the equivalent a 1 A — > 
AT 3 S~ ro-vibronic transition of O2 (u = p = 0) from 
Bellarv and Balasubramanianl (|1987f ). The individual 
line strengths are calculated for a given rotational tem- 



perature and then the profile is convolved to NIR- 
SPEC's resolution. We constructed a library of model 
spectra from 400K to 1500K at 50K intervals which 
we use to estimate the dominant rotational temper- 
ature of the emitting gas. In each case we use chi- 
squared minimization to obtain the best fit to the 
data, and from this we were able to estimate the er- 
rors shown in Fig. 01 We used data in the range of 
1.70 to 1.72 pm to fit the line, as this contains the 
region of highest signal to noise and largest variation 
with temperature. The results, shown superimposed 
in Fig. IH give rotational temperatures between 550K 
and 1000K, thus supporti ng the volcanic origin (DP02, 
Zolotov and Feglevl fl998h of the SO gas. The errors 
are of order 30-100K justifying our choice of 50K in- 
crements for our model spectra. 

At ~i.69/im there is also an unidentified peak 
which, as mentioned in DP02, corresponds to the S R 
branch of the SO band; however, it is much weaker in 
the laboratory spectra, possibly implying some non- 
thermal distribution of excited SO states. While we 
find this intriguing, analysis of distribution of states 
of SO gas in non-LTE conditions is beyond the scope 
of this paper. 

3.3. Total Photon Flux 

The total observed photon flux can be calculated 
by integrating over the observed emission band (from 
1.685-1.725 ^m), giving us values of F b s = 0.051, 
0.031, 0.022, 0.035, and 0.017 photons/s/cm 2 for the 
1999,2000, 2002, 2003 and 2002 AO nights respec- 
tively. For easier comparison all data have been ad- 
justed to an observer to Io (geocentric) distance of 
4.08AU. We can convert these to the total number of 
photons emitted by SO if we assume isotropic emis- 
sion of the excited atoms in a plume, N pe = F }, s 4:irc(j , 
where di is the geocentric distance to Io, giving 2.4 x 
10 27 , 1.5 x 10 27 , 1.0 x 10 27 , 1.6 x 10 27 and 0.8 x 10 27 
photons/sec respectively (see Tab. [2]). Note that these 
values are independent of surface emitting area or dis- 
tance to the observer. Various excitation mechanisms 
for SO are discussed in DP02, in an attempt to convert 
these values to a observed column density. However, 
ultimately more information is needed about the local 
physics. Simple thermodynamic equilibrium calcula- 
tions are rendered invalid by the rarefied nature of 
Io's atmosphere, where low collision rates would make 
it difficult to ensure a Boltzmann population of states. 
Note the observed peak at 1.69/im, which may also 
indicate a non-equilibrium distribution of states. 
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1.66 1.68 1.70 1.72 1.74 1.76 1.66 1.68 1.70 1.72 1.74 1.76 

Wavelength(yiim) Wavelength(/u.m) 

Fig. 3. — Raw lab spectra showing the various transition branches of the SO 1.7/um band (left) and convolved with 
NIRSPEC resolution (right) 



4. Discussion 

The study of this ro-vibronic transition gives us 
a valuable tool to investigate SO, a species thought 
to originate from volcanic out-gassing (DP02). The 
exact composition and physical conditions of Io's at- 
mosphere are still poorly constrained and while it is 
agreed that both volcanic processes and sublimation 
contribute to the atmosphere, the details of it's minor 
components are unclear. In this study we fit rotational 
temperatures ranging from 550K-1000K to the emis- 
sion band of SO gas. These high rotational tempera- 
tures suggest SO to be of a volcanic origin. A volcanic 
or high pressure vent origin is plausible, because inter- 
conversion of rotational and translational energy only 
requires tens of collisions in comparison to ~ 10 8 colli- 
sions for interconversion of electro nic (a 1 ) and transla- 



tiona l energy in the O2 molecule ([Boodaghians et al 



1983), so one would not expect the SO (a 1 A) electronic 
state to be in thermodynamic equilibrium with SO's 
rotational levels in the tenuous plume atmospheres on 
Io. 

An interesting parameter which remains uncon- 
strained is the total SO column abundance, which 
could be determined if the ro-vibronic band B l Y> — > 
A 3 £~ of SO at 0.97 /im would b e measured simulta - 
neously with the 1.707/xm band ( Setzer et a l. 1999). 
During the March 2003 run we therefore attempted 
to take spectra with NIRSPEC 1, a filter which cov- 
ers the range 0.947-1.121 fim. Unfortunately we could 
not detect the 0.97 /im transition, which is weak and 
dominated by the strong glare of reflected sunlight 
from Jupiter at these shorter wavelengths. The non- 



uniformity of the contribution from Jupiter in our slit 
also prevents estimates of an upper limit on the ratio 
of the 0.97/im and 1.7/mi band. 

In July 2000, we also observed with NIRSPEC in 
high resolution echelle spectroscopy mode (R=25000), 
in an attempt to resolve individual transitions, but the 
1.7 /im band was not detectable. We conclude that 
a combination of limited observing time during Io's 
eclipse, the higher dispersion of our signal and a low 
band-to-noise ratio as seen in the low resolution data 
of the same time period, result in a signal to noise too 
low for us to detect. 

We are able to add to the estimates of emitted flux 
from the a 1 A — > A 3 £~ band, previously reported in 
DP02. The strongest emission is found in the 1999 
data, which also happens to be the observation dur- 
ing which Loki was b r ightes t, as shown in Fig. 1 of 
Rathbun and Spencer (2006). In Table [2] we summa- 



rize the SO emission, rotational temperature and the 
3.5/im flux measured from Loki (J. Rathbun, personal 
communication). The highest photon flux is found in 
the 1999 data, which also happens to be the observa- 
tion when Loki was brightest. Although naively one 
might expect the high photon fluxes to also correspond 
to high temperatures, no such correlation is seen. The 
temperature in 1999 is high, 900±30K, but a similar 
high temperature is seen in Nov 2002 when the photon 
flux was lowest. This may be caused by the fact that 
we may observe several active volcanic centers simulta- 
neously, and our non-AO observations cannot resolve 
individual volcanoes. In addition, usually only a few 
of the hot spots are associated with plumes and the 
temperature in these plumes can vary from volcano to 
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Fig. 4. — SO band emission for all nights. Over-plotted are the best fit rotational temperature models using chi- 
squared minimization. The errors quoted in the legend the 1 sigma uncertainties. Also shown are example line 
profiles at 322K and 1500K for comparison. 
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Table 2 

Integrated band emission of SO's a 1 A — * A 3 £~ transition from 1999-2003 



Date 


SO emission 


Rotational Temp. 


Loki Measurement Date b 


Loki Flux b 




(xlO 27 Photons/sec) 


(K) 




(GW/^m/str) 


9/24/1999 


2.4 ± 0.7 


900 ± 30 


9/15/1999 - 10/10/1999 


22 - 60 


7/19/2000 


1.5 ± 0.5 


550 ± 80 






11/12/2002 


0.8 ± 0.2 a 


850 ± 90 


10/27/2002 


6 


12/21/2002 


1.0 ± 0.3 


550 ± 100 


10/27/2002 - 1/13/2003 


6 - 29 


3/25/2003 


1.6 ± 0.5 


1000 ± 70 


3/25/2003 


30 



a Disk averaged AO data 

b Data provided by Julie Rathbun (private communication). In each case we have shown the Loki observations 
that are closest in time to our SO data. No 3.5/Ltm measurements of Loki were taken within 4 months on either side 
of the 19th July 2000 observations. 



volcano. We see the average effect of all of these. 

With so few data points it is difficult to draw any 
solid conclusions. However, it should be noted that the 
highest band emission occured in late 1999, concur- 
rent wit h a dramatic rise in Loki's o bserved 3.5/zm in- 
tensity (jRathbun and Spencer 2006). Also the lowest 
observed SO band emission of 0.8 x 10 27 photons/sec 
occurs in late 2002 which corresponds to a time of very 
low Loki flux. 

Accurate measurements of the background IR con- 
tinuum would allow us to determine the behavior of the 
SO emission in the context of the broader thermal IR 
emission of Io. However further study of the raw spec- 
tra indicates that contamination from reflected and 
diffuse Jupiter light affect the background level and 
slope, particularly in the 2000 and 2002 data, due to 
observational factors such as slit orientation and prox- 
imity to Jupiter. By their nature these observations 
occur when Io is very close to Jupiter (20-40"), and 
that distance is changing rapidly in time resulting in 
a varying level of contamination from frame to frame. 
It should be noted that this does not affect our mea- 
surements of the SO emission band. 



Comparison of the band strengths on each of the 
nights show some variation. Although the photo- 
metric calibration is rough there appears to be a 
genuine change in SO emission from the 1999 data 
(2.4 x 10 27 photons/sec) to subsequent nights (~1 x 
10 27 photons/sec). This is mirrored by the measured 
flux of Loki over this time period indicating that SO 
emission may be linked with Loki's activity, although 
further investigation at high spatial resolution and/or 
concurrent monitoring of Loki and disk integrated SO 
spectra is required to prove there is a definitive corre- 
lation. 

Io's spectrum, when the satellite is in eclipse, is de- 
termined by hotspots, and so we expect to see a large 
variation in band-to-continuum ratio of SO as we look 
at finer spatial scales on Io. Using AO techniques one 
can spatially resolve individual hotspots, and perhaps 
be able to confirm the confinement of the SO 1.707 um 
emission band to volcanic areas. Ultimately, we hope 
to determine the exact spatial distribution of SO in Io's 
atmosphere, to determine whether Io's atmosphere is 
supported predominantly by volcanic activity or sub- 
limation of SO2 frost. 



5. Conclusion 

Sulfur monoxide was detected in the near infrared 
at 1.707 fan in all of the low resolution NIRSPEC 6 
data sets taken from 2000 to 2003. Analysis of the 
NIRSPEC 1 data taken on 25th March 2003 failed 
to detect the 0.97 fim band, ( Setzer et al.lll999l h due 
to contamination by reflected sunlight from Jupiter. 
High spectral resolution data taken in July 2000, also 
failed to resolve the 1.7/im band. 



6. Acknowledgments 

We would like to extend our warmest thanks to 
Julie Rathbun, for providing the data on Loki's ac- 
tivity. We also thank Emmanuel Lellouch and John 
Pearl for their very constructive and helpful reviews. 
Support for this research was provided by NSF grant 
AST 0406275 and by the HCICT Program (NASA 
Grant No. NNG05GA25G). The data were obtained 
with the W.M. Keck Observatory, which is operated 
by the California Institute of Technology, the Univer- 



8 



sity of California, Berkeley and the National Aero- 
nautics and Space Administration. The Observatory 
was made possible by the generous financial support 
of the W.M. Keck Foundation. The authors extend 
special thanks to those of Hawaiian ancestry on whose 
sacred mountain we are privileged to be guests. With- 
out their generous hospitality, none of the observations 
presented would have been possible. 

REFERENCES 

Ballester, G. E., McGrath, M. A., Stobel, D. F., Zhu, 
X., Fcldman, P. D., and Moos, H. W. (1994). Detec- 
tion of the S02 atmosphere on Io with the Hubble 
Space Telescope. Icarus, 111:2-17. 

Bellary, V. P. and Balasubramanian, T. K. (1987). 
On the rotational intensity distribution in the 
a 1 A g — V 3 £~ magnetic dipole transition of oxy- 
gen molecule. Journal of Molecular Spectroscopy, 
126:436-442. 

Boodaghians, R. B., Borrell, P. M., and Bor- 
rell, P. (1983). Room-temperature rate con- 
stants for the gas-phase quenching of metastable 
molecular oxygen, 2 (a 1 A g ) and 2 (b 1 'E g ~), by 
C0 2 ,N 2 0,NO,NH 3 ,HCl and S0 2 . Chemical 
Physics Letters, 97:193-197. 

de Pater, I., Laver, O, Marchis, F., Roe, H. C, and 
Macintosh, B. A. (2007). Spatially Resolved Ob- 
servations of the Forbidden SO a 1 A - A 3 £~ Rovi- 
bronic Transition on Io during an Eclipse . Icarus 
(Submitted). 

de Pater, I., Marchis, F., Macintosh, B. A., Roc, H. C, 
Le Mignant, D., Graham, J. R., and Davies, A. G. 
(2004). Keck AO observations of Io in and out of 
eclipse. Icarus, 169:250-263. 



Lellouch, E., Strobel, D. F., Belton, M. J. S., Sum- 
mers, M. E., Paubert, G., and Moreno, R. (1996). 
Detection of Sulfur Monoxide in Io's Atmosphere. 
ApJ, 459T107-110. 

Macintosh, B. A., Gavel, D., Gibbard, S. G., Max, 
C. E., Eckart, M., dc Pater, I., Ghez, A. M., and 
Spencer, J. (2003). Speckle imaging of volcanic 
hotspots on Io with the Keck telescope. Icarus, 
165:137-143. 

Marchis, F., de Pater, I., Davies, A. G., Roe, H. G., 
Fusco, T., Lc Mignant, D., Descamps, P., Macin- 
tosh, B. A., and Prange, R. (2002). High-Resolution 
Keck Adaptive Optics Imaging of Violent Volcanic 
Activity on Io. Icarus, 160:124-131. 

Marchis, F., Le Mignant, D., Chaffee, F. H., Davies, 
A. G., Kwok, S. H., Prange, R., de Pater, I., Amico, 
P., Campbell, R., Fusco, T., Goodrich, R. W., and 
Conrad, A. (2005). Keck AO survey of Io global vol- 
canic activity between 2 and 5 ^m. Icarus, 176:96- 
122. 

McGrath, M. A., Belton, M. J. S., Spencer, J. R., 
and Sartoretti, P. (2000). Spatially Resolved Spec- 
troscopy of Io's Pele Plume and S0 2 Atmosphere. 
Icarus, 146:476-493. 

Pearl, J., Hanel, R., Kunde, V., Maguire, W., Fox, 
K., Gupta, S., Ponnampcruma, C, and Raulin, F. 
(1979). Identification of gaseous S02 and new upper 
limits for other gases on Io. Nature, 280:755-758. 

Persson, S. E., Murphy, D. C, Krzeminski, W., Roth, 
M., and Rieke, M. J. (1998). A New System of Faint 
Near-Infrared Standard Stars. The Astronomical 
Journal, 116:2475-2488. 



de Pater, I., Roe, H., Graham, J. R., Strobel, D. F., 
and Bernath, P. (2002). NOTE: Detection of the 
Forbidden SO a 1 A - V 3 XV Rovibronic Transition 
on Io at 1.7 fim. Icarus, 156:296-301. 

Ingcrsoll, A. P. (1989). Io meteorology - How atmo- 
spheric pressure is controlled locally by volcanos 
and surface frosts. Icarus, 81:298-313. 

Kliore, A., Cain, D. L., Fjeldbo, G., Seidel, B. L., and 
Rasool, S. I. (1974). Preliminary Results on the 
Atmospheres of Io and Jupiter from the Pioneer 10 
S-Band Occultation Experiment. Science, 183:323- 
324. 

Lellouch, E., Encrcnaz, T., Belton, M., de Pater, I., 
and Gulkis, S. (1990). Io's atmosphere from mi- 
crowave detection S02. Nature, 346:639-641. 



Rathbun, J. A. and Spencer, J. R. (2006). Loki, Io: 
New ground-based observations and a model de- 
scribing the change from periodic overturn. Geo- 
phys. Res. Lett, 33:17201. 

Setzer, K. D., Fink, E. H., and Ramsay, D. A. (1999). 
High-Resolution Fourier- Transform Study of the 
- V 3 £- and a 1 A - X 3 T,~ Transitions of SO. 
Journal of Molecular Spectroscopy, 198:163-174. 

Strobel, D. F., Zhu, X., and Summers, M. F. (1994). 
On the vertical thermal structure of Io's atmo- 
sphere. Icarus, 111:18-30. 

Summers, M. E. and Strobel, D. F. (1996). Pho- 
tochemistry and Vertical Transport in Io's Atmo- 
sphere and Ionosphere. Icarus, 120:290-316. 



9 



Wong, M. C. and Johnson, R. E. (1996). A three- 
dimensional azimuthally symmetric model atmo- 
sphere for Io 1. Photochemistry and the accumu- 
lation of a nightside atmosphere. J. Geophys. Res., 
101:23243-23060. 

Wong, M. C. and Smyth, W. H. (2000). Model Calcu- 
lations for Io's Atmosphere at Eastern and Western 
Elongations. Icarus, 146:60-74. 

Zhang, J., Goldstein, D. B., Varghese, P. L., 
Gimelshcin, N. E., Gimclshcin, S. F., and Levin, 
D. A. (2003). Simulation of gas dynamics and radi- 
ation in volcanic plumes on Io. Icarus, 163:182-197. 

Zolotov, M. Y. and Fegley, B. (1998). NOTE: Volcanic 
Production of Sulfur Monoxide (SO) on Io. Icarus, 
132:431-434. 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.2. 



10 



